home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d18 / nrpas13.arc / FTEST.PAS < prev    next >
Pascal/Delphi Source File  |  1991-05-01  |  915b  |  31 lines
  1. PROCEDURE ftest(data1: gln1array; n1: integer;
  2.       data2: gln2array; n2: integer; VAR f,prob: real);
  3. (* Programs using routine FTEST must define the types
  4. TYPE
  5.    gln1array = ARRAY [1..n1] OF real;
  6.    gln2array = ARRAY [1..n2] OF real;
  7. in the main routine, with n2 less than or equal to n1. *)
  8. VAR
  9.    i: integer;
  10.    var2,var1,df2,df1,ave2,ave1: real;
  11. BEGIN
  12.    IF (n2 > n1) THEN BEGIN
  13.       writeln('pause in routine FTEST');
  14.       writeln('first array must be the larger'); readln
  15.    END;
  16.    avevar(data1,n1,ave1,var1);
  17.    FOR i := 1 TO n2 DO data1[i] := data2[i];
  18.    avevar(data1,n2,ave2,var2);
  19.    IF (var1 > var2) THEN BEGIN
  20.       f := var1/var2;
  21.       df1 := n1-1;
  22.       df2 := n2-1
  23.    END ELSE BEGIN
  24.       f := var2/var1;
  25.       df1 := n2-1;
  26.       df2 := n1-1
  27.    END;
  28.    prob  :=  betai(0.5*df2,0.5*df1,df2/(df2+df1*f))
  29.          +(1.0-betai(0.5*df1,0.5*df2,df1/(df1+df2/f)))
  30. END;
  31.